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The interference pattern of the resonance fluorescence from aJ = l/2toJ = l/2 transition 
of two identical atoms confined in a three-dimensional harmonic potential is calculated. Thermal 
motion of the atoms is included. Agreement is obtained with experiments [Eichmann et ai, Phys. 
Rev. Lett. 70, 2359 (1993)]. Contrary to some theoretical predictions, but in agreement with the 
present calculations, a fringe visibility greater than 50% can be observed with polarization-selective 
detection. The dependence of the fringe visibility on polarization has a simple interpretation, based 
on whether or not it is possible in principle to determine which atom emitted the photon. 
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I. INTRODUCTION 

Many variants of two-slit interference experiments, often 
"thought experiments," have been used to illustrate funda- 
mental principles of quantum mechanics. Recently, Eichmann 
et al. ||lj have observed interference fringes in the resonance 
fluorescence of two trapped ions, analogous to those seen in 
Young's two-slit experiment. Of particular interest was the 
fact that the interference fringes appeared when it was im- 
possible in principle to determine which ion which scattered 
the photon disappeared when it was possible. This is in 
agreement with Bohr's principle of complementarity, which 
requires that the wave nature of the photon (the interference 
fringes) cannot be observed under the same conditions as its 
particle nature (the possibility of assigning to the photon a 
trajectory that intersects just one of the ions). In contrast 
to many thought experiments the disappearance of the 
fringes when the path of the particle can be determined has 
nothing to do with the position-momentum indeterminacy re- 
lations. The experiment contains features from some thought 
experiments of Scully and Driihl j^, regarding the interfer- 
ence of light scattered by two multi-level atoms. 

Recently, controversy has arisen over the mechanism by 
which complementarity is enforced in a two-slit interference 
experiment. Some claim that the destruction of interference 
by a determination of the particle's path is always due to a 
random momentum transfer necessitated by the indetermi- 
nacy relations Others claim that the mere existence 
of the path information can be sufBcient to destroy the in- 
terference [^]. Englert et al. claim that the experiment of 
Eichmann et al. supports the second position 1^. 

Published calculations explain some aspects of the obser- 
vations of Eichmann et al. |p|-[l^. However, none of those 
calculations include all of the factors required to make a com- 
parison with the experimental data. Here, we calculate the 
scattering cross section for arbitrary directions and polariza- 
tions of the incident and outgoing light. While the results 



were used in the analysis of the data in Ref. |l|, the details 
of the calculations were not given. The main limitation of 
the calculation is the use of perturbation theory, so that it is 
valid only for low light intensities. However, it includes the 
effect of thermal motion more precisely than has been done 
elsewhere, taking into account the actual normal modes of 
the system. Also, the actual experimental geometry is fully 
taken into account, which is not always the case in the other 
calculations. 

Finally, we clarify the sense in which the loss of the fringe 
visibility [defined as (J^ax - /min)/(/max + /min)] for certain 
detected polarizations is due to the existence of "which path" 
information in the ions. This is an application of the funda- 
mental quantum principle that transition amplitudes are to 
be added before squaring if and only if they connect the same 
initial and final states. 



II. EXPERIMENT 

The experimental apparatus has been described previously 
. Figure ^ shows the geometry. Two ^^®Hg^ ions were 
confined in a linear Paul (rf) trap by a combination of static 
and rf electric fields. The ions were laser-cooled to tem- 
peratures of a few mK with a beam of linearly polarized, 
continuous- wave light, nearly resonant with the 194 nm tran- 
sition from the ground 6s ^51/2 level to the 6p^Pi/2 level. The 
laser beam diameter was about 50 fim, and the power was 50 
fiW or less. The same beam was the coherent source for the 
Young's interference. Cooling in the trap resulted in strong 
localization of the ions, which was essential for observation 
of interference fringes. The trap potentials were arranged so 
that a pair of ions would be oriented along the symmetry 
(Z) axis of the trap. The incoming photons, with wavevec- 
tor kin and polarization vector ein, made an angle Q of 62° 
with respect to the Z axis. The X axis is oriented so that 
the X-Z plane contains kin. Light emitted by the ions was 
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coUimated by a lens and directed to the surface of an imaging 
photodetector, which was used to observe the fringes. The 
wavevector and polarization of an outgoing photon are kout 
and eout- The projection of kout onto the X-Z plane makes 
an angle with respect to kin. The deviation of kout from 
the X-Z plane in the +y direction is $ (not shown in Fig. |l|) . 
The sensitive area of the photodetector included a range of 
i\> from about 15° to 45° and a range of $ from about -15° 
to +15°. For polarization-selective detection, a glass plate 
oriented at Brewster's angle was placed in the detection path 
so that nearly all of the light with eout in the X-Z plane was 
transmitted into the glass, while some of the light polarized 
along the Y axis was reflected to the imaging detector. The 
input polarization ein was varied. Another lens system formed 
a real image of the ions on a second imaging detector. This 
image was used to determine when there were precisely two 
ions in the trap. 



The number operators are defined in the usual way by Nf°'^ = 
(a°°™)ta°°"" and A^f' = (a^°')ta^°i for i = X,Y,Z. The anni- 
hilation operators are defined in the usual way, for example: 

com / TtlUJz com . ^ 7-jcom / p- \ 

az =\/^— + Pz ■ (5) 

The three center-of-mass modes have the same frequencies as 
those of a single ion in the trap: ljz and ljh. The three relative 
modes include a symmetric stretch mode along the Z direction 
at frequency los = \Pi<~oz and two tilting or rocking modes 
along the X and Y directions at frequency ujt = (uJ\i — lo\')^^'^ 
The eigenstates of -fftrans are the simultaneous eigenstates of 
the set of number operators | n'^"", n^°"", nf"", n5?', n?=', nf ') 
with eigenvalues ft[wz(nf"" + 1/2) + ^^RinT" + n-y™ + 1) + 
us(n'i 1/2) + ujTin'i + n^°' + 1)]. 



IV. ATOMIC LEVEL STRUCTURE 



III. TWO-ION HARMONIC OSCILLATOR 
SYSTEM 



In the pseudo-potential approximation, the Hamiltonian for 
the translational motion of the two ions in the harmonic trap 
is 

Htrans = ^ + ^ + V"(Rl) + V(R2) + 



2m 2m 



47reo I Ri - R2 



(1) 



where Ri and Pi are the position and momentum of the ith 
ion, e and m are the charge and mass of an ion, and 



1/(R) = \nu4i {X" + F") + \nuJ-zZ'' 



(2) 



is the potential energy of a single ion in the trap. In Eq. (Q), 
we have made the approximation that the trap pseudopo- 
tential is cylindrically symmetric. Here, R=(X, F, Z), in the 
Cartesian coordinate system shown in Fig. ly. The classical 
equilibrium positions of the ions, found by minimizing the to- 
tal potential energy, are R? = (d/2)Z and R" = -(d/2)Z, 
where d — [e^ /{2-K€(,mijLi%)Y^^ , and it is assumed that ljr > 

Ulz- 

For small displacements ui — Ri — R^ and U2 = R2 — R2 
about the equilibrium positions, a harmonic approximation 
can be made. The Hamiltonian separates into terms involv- 
ing either center-of-mass (com) or relative (rel) coordinates 
and momenta defined by 



U°°"^ = (Ui + U2)/2 

u''°' = (Ui - U2)/2 

pcom = ^ 

p-1 = _ P2 



(3) 



The translational Hamiltonian, in the harmonic approxima- 
tion, is 

Htrans = hiOziN^z'^ + 1/2) + hLOR{N^^^ + + 1) 

+ huJs{Nt + 1/2) + fu^T{Nt + iVy ' + 1). (4) 



Figure g shows the magnetic sublevels involved in the 
6s ^51/2 to 6p^Pi/2 transition. These levels form an ap- 
proximately closed system, since the probability that the 
&p^P\/2 level radiatively decays to the 5d^6s^ ^0^/2 level is 
only 1.4 X 10~^ |l^. The rest of the time it returns to the 
ground 6s ^51/2 level. The 5d®6s^ ^D2,/2 level has a lifetime of 
9 ms and decays with about equal probability to the ground 
level or to the 5ci^6s^ ^^5/2) which has a lifetime of 86 ms and 
decays only to the ground level. 

Since the static magnetic field is small, we are free to define 
the quantization axis of the ions to be along the electric polar- 
ization vector ein of the incident light. If the static magnetic 
field is along some other direction, then the Zeeman sublevels 
defined according to the electric polarization vector are not 
stationary states. This does not change the analysis as long 
as the Zeeman precession frequency is much less than the in- 
verse of the scattering time, which is approximately equal to 
the 6p ^Pi/2-state lifetime (2.3 ns). In the experiments de- 
scribed here, the magnetic field was small enough that this 
was always the case. 

Figure ^ shows a Cartesian coordinate system having its 
z axis oriented along ein. The x axis is parallel to kin. The 
y axis is defined so that (x, y, z) forms a right-handed coor- 
dinate system. This coordinate system is more useful than 
the trap-oriented (X, Y, Z) coordinate system of Fig. |l] for 
describing the angular distribution of the scattered light. 



V. SCATTERING CROSS SECTION 

Consider the process in which two ions, initially in their 
ground electronic states, absorb a photon having wavevector 
kin and polarization ein, emit a photon having wavevector 
kout and polarization eout, and are left in their ground elec- 
tronic states. The ions may change their Zeeman sublevels 
during the process. Also, the motional state of the two ions 
may change. 

The electric dipole Hamiltonian which causes the transi- 
tions is 



7?ED = -Di • E(Ri,t) - D2 • E(R2,t), 



(6) 
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where Di and D2 are the electric dipole moment operators 
for ions 1 and 2 and E(R, t) is the electric field, consisting of 
a classical part, representing the incident laser beam, and the 
quantized free field operator: 



hujs 
2^ 



(7) 



where 5R denotes the real part, fo is the amplitude of the laser 
electric field, and as is the annihilation operator for a photon 
of wavevector ks, frequency ujs, and polarization Ss, and V is 
the quantization volume. 

The electric dipole Hamiltonian, in second-order perturba- 
tion theory, gives the cross section for the two ions to scatter 
a photon in a particular direction: 



|*,){*,|(Di.ei„)e*--'^ilvl>,) 



ujo - t^iii + {Ej - Ei)/h- Z7/2 



UJO - i^in + {Ej - Ei)/Ti - 17/2 



(8) 



where LO'ia — cj kin 

I, ojout ~ c| kout I, TiLoo is the separa- 
tion between the ground and excited electronic states of an 
ion, 7 is the decay rate of the excited state, and Ci = 
'^out/(167r^c''?i^eo)- The initial, final, and intermediate states 
describing the electronic and motional degrees of freedom of 
the system are l^^i), 1*1//), and The energies Ei, Ef, and 

Ej are the motional energies of the ions in the initial, final, 
and intermediate states. They depend on the values of the 
six harmonic oscillator quantum numbers, which we denote 
by {jiho}. Because of energy conservation, the frequency of 
the outgoing photon depends on the final state: 

t.^'out = LOin + (Ei ~ Ef)/h. (9) 

Thus, the scattered light has a discrete frequency spectrum, 
and the different components could, in principle, be detected 
separately. In Eq. (^, all frequency components are summed, 
which is appropriate if the detection is frequency-insensitive. 
The laser frequency is assumed to be nearly resonant with an 
optical transition in the ion, so that only one intermediate 
electronic state has to be included in the sums, and we can 
neglect the counter-rotating terms. We ignore dipole-dipole 
interactions between the ions, because they were separated 
by many wavelengths in the experiment. Here we specialize 
to the case of an ion which has no nuclear spin and which 
has a ^51/2 ground state and a excited state, like the 

i^^Hg"*" ions used in Ref. Q. We denote a state in which ion 
1 is in the (^51/2, mj = -1-1/2) state, ion 2 is in the (^Pi/2, 
mj — —1/2) state, and has the harmonic oscillator quantum 
numbers {riho} by 

1*) = |('Si/2, +1/2)1(^^/2, -l/2)2{nho}). (10) 

There are four possible sets of initial mj quantum numbers 
for the two ions and four possible final sets. There are two ba- 
sic kinds of scattering processes — those which preserve the 
nij quantum numbers of the ions, and those which change 
mj of one ion. We treat these cases separately. The form of 
Eq. (^ excludes the possibility of both ions changing their 
mj quantum numbers. 



A. Both mj quantum numbers remain the same (tt 
case) 

In order to be definite, we let mj — +1/2 for both ions, 
both before and after the scattering. That is, 

I*.) = i('Si/2, +1/2)1(^51/2, +l/2)2{nhoW, (11) 
1*/) = l('Si/2,+l/2)i(^Si/2,+l/2)2{nho}/). (12) 

We call this the vr case, because it involves only n transitions, 
that is, transitions that leave mj unchanged. Because of the 
electric dipole selection rules, the only intermediate states 
which contribute nonzero terms are of the form 

1*.) = l('Pi/2,+l/2)i('Si/2,+l/2)2{nho},) (13) 

for the first sum, and 

I*.) = l('Si/2,+l/2)i(2Pi/2,+l/2)2{nho},) (14) 

for the second sum. The matrix elements connecting the ini- 
tial states to the intermediate states are 

{*,i(D,.ei„)e*--'^<'|*,) 
= {('Pi/2,+l/2)piDj,,!(^Si/2,+l/2)p) 

x{{nho},|e*'"-^''i{nho}>) 
= -l=(^Pi/2||D(^'f Si/2){{nhob|e*-^p|{nho}.), (15) 

where p = 1 or 2, Dpz is the z component of the Dp operator, 
and (^Pi/2 ll^'^' lP'S'1/2) is the reduced matrix element of the 
dipole moment operator (the same for both ions). 

The angular distribution of the outgoing photon is con- 
tained in the matrix elements connecting the intermediate 
states to the final states. The unit propagation vector for the 
outgoing photon is 

kout ~ (sini? cos sin I? sin (p, cos I?), (16) 

where ■& and ip are spherical polar angles with respect to the 
{x,y,z) coordinate system of Fig. ^. The polarization vector 
lout must be perpendicular to kout. We define two mutually 
orthogonal unit polarization vectors, both perpendicular to 
kout, by 
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and 



— {— cos 1? COS If, — COS i9 sin Lp, sin i?) 



= (— sin ((9, cos ((9, 0). 



(17) 



(18) 



Since only the z components of Di and D2 contribute to the 
matrix elements connecting the intermediate and final states 
for this case, light with polarization vector icr cannot be emit- 
ted. 



With the choice of eout=e7r, the matrix elements connecting 
the intermediate states to the final states are 

{*/|(Dp.e.)e-*°--^-!*,) 
= sin^{('Si/2,+l/2)j,|Z)p,!('Pi/2,+l/2)p) 



x{{nh4^ie-*--^-|{nhob) 



sin i9 



(^5i/2||D(i'fPi/2){{nho}/|e- 



-ikout -R. 



-\{n^o},). (19) 



Equation (jsj) for the cross section becomes 
dcr*^' sin^ ■& , 



36 

+ E 

i 



E 



{{nho}/|e 



|{nho},){{nho},|e"'-^M{^ho}.) 



{{nho}/|e 



-ikout-R-s 



|{nho}j){{nho}i|e" 



ojQ - Ci;in + {Ej - Ei)/h- 17/2 



H{nho}i) 



LOQ - tJin + (-Bj - Ei)/h - 17/2 



(20) 



The (1) superscript on the cross section is to label it as the 
TT case. The same result would have been obtained for any of 
the other three possible sets of initial mj quantum numbers, 
so this is the general result for the case in which the m,j values 
do not change. The presence of two terms in Eq. ( |2o| , which 
are added and then squared, is the source of the Young's in- 
terference fringes. These two terms can be identified with the 
two possible paths for the photon, each intersecting one of the 
two ions. 

The sums over intermediate harmonic oscillator states can 
be done by closure if the energy denominators are constant. 
While they are not constant, because [Ej — Ei) varies, it can 



be shown (see, for example, Ref. \^) that the main contri- 
butions to the sum come from terms where \Ej — Ei\ is less 
than or on the order of ^ REi, where R is the photon re- 
coil energy (?ifcout)^/(2m). For the Hg^ 194.2 nm transition, 
R—h X 26.7 kHz. For Doppler cooling, Ei will be on the order 
of hy, where, for this transition, h'y—h x 70 MHz. Thus, the 
rms value of (Ej — Ei) will be on the order of ^ x 1.4 MHz, 
which is much less than h'y/2. Therefore, while the denomi- 
nators are not strictly constant, they are nearly constant for 
the terms which contribute significantly to the sums. 

If we neglect {Ej — Ei)/h compared to 7/2 and use closure 
to evaluate the sums, Eq. (fed) simplifies to 



dcr^^^ sin^ 1? 1(^51/211-0^^^ ||^Pi/2)l'' v;^ 1,, , , -ikoutRi iki„.R, , r 1 \ , /r 11 -^out-Ro iki„.R,ir 1 \| 
dilout 36 (ljo — t^in) + 7 /4 ^ ' ' 

Sin^l9 !(^Si/2||-D*^'f Pl/2)|* „ 1,, , , -iq.Rnr 1 \ , /r l 1 -iq-R2ir 1 \|2 

= -^^ (cjo -a;in)2 +774 Z^'^M<{"-''°^-f!^ MWo}») + {{nho}/|e 'l{nho}>)| 



'^|(^Si/2||I>«fPl 



/2i 



36 (ljo - t^in)2 7V4 



(21) 



where q = kout — kin. Since the branching ratio for decay of 
the excited '^Pi/2 states to the ground states is nearly 

100%, the spontaneous decay rate 7 is 



7 = 



Equation ( |2l| ) for the cross section becomes 



da^^^ _ sin^i9 

xE 



(22) 



(23) 



where ctq = Ao/27r is the resonance cross section, Ao = 2-kc/loo 
is the resonance wavelength, and £(cjin — 1^0) is a Lorentzian 
of unit height and width 7: 



(7/2) = 



(o-in - UJoY + (7/2)2 ■ 



(24) 



In deriving Eq. (|2^), we have assumed that dt;o/i.^out ~ 1. The 
sum over final harmonic oscillator states can be done by clo- 
sure: 



dcrt^> _ sin^-i? 

_ sin^ d 
~ Stt 



aoC{uJin — t^o)({who}il(( 



iqR]_ . „iqR: 



+ e" 



+ e 



-iq-R,2 



)j{nho}i) 



iq.(Rl-R2) I -iq.{Ri-R2) 



(25) 
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The exponentials can be combined in Eq. ( pq ) because the components of Ri and R2 commute. The cross section can be 
written in terms of the equilibrium ion separation d and the displacement coordinates ui and U2 as 



dcr^i' sin^ i9 



Stt 



ao£(a;i„-c^o)({nho},l2 + e"''-('*+"i-"^)+e-''>''<*+"i-"^'|{nho},:)- 



(26) 



The exponential factors in Eq. ( |26[ ) depend on the relative coordinates of the two ions and not on their center-of-mass coordi- 
nates. 

In order to compare with the experiment, we compute the cross section averaged over a thermal distribution of |{nho}i) 
initial states: 



Stt 



-cro£{uJin — uJo) [2 + e''^ '^ (e 



iq'(ui-U2) 



(27) 



where (A) denotes the thermal average of the operator A. For harmonic oscillators, the thermal averages have a simple form 

TIF 



(e 



±iq-(ui-U2) 



i([q.(ui-U2)]^> 



(28) 



While Refs. assume a common temperature for all of the harmonic oscillator modes, Eq. ( p8| ) is still valid if different 

modes have different temperatures. Different modes are laser-cooled at different rates depending on the direction of the laser 
beam. Hence, the modes have different temperatures unless the energy transfer rate between them is fast [Q. The thermally 
averaged cross section is 



da 



(1) 



47r 



aoC{LL>in — uJo) l + cos(q-d)e 2 



i([q-(ui-U2)]2> 



(29) 



which is equivalent to Eq. (I) of Ref. except that it includes the sin i9 angular dependence. The interference fringe visibility 
is given by the exponential factor multiplying cos(q • d) This factor decreases with increasing temperature and is analogous to 
the Debye- Waller factor for x ray scattering from a crystal. It can be rewritten as 



4([q-(>^ 



2)] = > 



: exp 



■ exp 



: exp 



Wx 
2mujT 



\ 2/ mujT \ 2/ mijjs \ 2/ 



coth 



l^2fcflrf=i 



2mujT 



coth 



mtos 

I hu)T 

\2kBT{ 



rrrcl 



2mu)s 



coth 



(30) 



where r|-° is the temperature of the mode, etc., and the 
approximation in the last line is valid when the mean har- 
monic oscillator quantum numbers are large. In the limit of 
small thermal motion or small |q| (near- forward scattering), 
the visibility can approach 100% (with polarized detection), 
in agreement with Ref. but in contradiction to Ref. 
where it was claimed that the visibility could not exceed 50%. 

B. One mj quantum number changes (cr case) 

Here we consider the case in which one of the ions changes 
its m,j quantum number in the scattering process. We call 
this the a case, since it involves a a transition, that is, a tran- 
sition that changes raj by ±1 in one of the ions. There are 
eight cases, since there are four possible initial states and two 
ions which could change quantum numbers. 

In order to be definite, we pick the case where mj = -1-1/2 
for both ions before the scattering, and ion 1 changes to 
mj = —1/2 after the scattering. That is, 

1*0 = i('5i/2,+l/2)i(2Si/2,+l/2)2{nho}.), (31) 



1*/) = i('5i/2,-l/2)i(^Si/2,+l/2)2{nho}/). (32) 

Only the first sum over j in Eq. (^) contributes, since only it 
contains Di, the dipole moment which leads to the change in 
mj of ion 1. 

As in the previous case, the only intermediate states which 
contribute nonzero terms are of the form 

= i('j'i/2,+l/2)i('Si/2,+l/2)2{nho},). (33) 

The matrix elements connecting the initial states to the in- 
termediate states are 

{vl/,j(Di-ein)e*--^i|*,) 
= (('Pi/2,+l/2)i|i?i.|('Si/2, +1/2)1) 

= -^(^Pi/2|l^<'^||'5i/2){{rzhob|e*--^MWo}.). (34) 

In the i9 = 7r/2 plane, only the polarization corresponding 
to Co- is emitted, but in general, light with both and ejr 
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contributes to the scattered intensity. We consider these two 
cases separately. 

For eout=e<T, the matrix elements connecting the interme- 
diate states to the final states are 



{*/|(Di.e%)e- 



-{(^Si/2,-l/2)i|i3(L\l(^Pi/2,+l/2)i) 



V2 



x({nho}/i£ 



V6 



-Cs 



1/2 



\D 



(1)||2 



{nho}j) 
-Pi/2)({nho}/ 



(35) 



where D^}^ is the (1,-1) spherical tensor component of the 
dipole moment operator for ion p. The rest of the calculation 
is very similar to the it case. The final result, analogous to 
Eq. ^2% for the n case, is 



Stt 



(36) 



which is independent of kout and shows no interference fringes. 
The same result would have been obtained for any of the other 
three initial states, since the absolute squares of the matrix 
elements are the same. 

For eout=e7r, the matrix elements connecting the interme- 
diate states to the final states are 



(vI>/|(Di-e,)e 

-cos'de''^ 



-ikout-R-l j^jy ,^ 



V2 

x({"-ho}/|e 
-cosiJe''^ ,2 



((^Si/2,-l/2)ilDii\l(^Pi/2,+l/2)i) 
""""'■''Mlnho},) (37) 



CS^M\D''''fPi/2){{n^o}f\^ 



V6 

The final result is 



cos^ ■& 
Sir 



aoiC{uJin — too), 



^'\{nho}j)- 



(38) 



which shows no interference fringes. The same result would 
have been obtained for any of the other initial states. 

For the a case, Young's interference fringes are not ob- 
served because only one of the two terms inside the absolute 
value bars in Eq. (H) is nonzero. There is only one path for 
the photon, intersecting the ion whose state is changed in the 
scattering process. 



C. Total cross section with or without 
polarization-selective detection 

In Ref. [Ql, a linear polarizer was sometimes placed before 
the photon detector. For experimental convenience, the orien- 
tation of this polarizer was fixed, while the input polarization 
could be varied. To obtain the total cross section describing 
a given experimental situation, we sum over all final atomic 
states and average over all initial states. For polarization- 
insensitive detection, we also sum over the polarizations of 
the outgoing photon. 

The cross section for polarization-insensitive detection is 



da 



unpol 



da 



(1) 



/ da(3) 



ao_ 

4% 



+ 2{ — ) + 2 ( 

C(uJin — CJo) { 1 + COS^ l9 

l + cos(q-d)e-^<['i'("i-"2))') 



sin'^ j9 



}• (39) 



The fringe visibility in this case cannot exceed 50%. 

The cross section for detection of light with polarization e.^ 

is 



da(i) \ ^ / da^'> 
+ 2 



da^ 

driout I \ d^out / \ di^out 
= ^L(y)\i^ - i^o) {cos^ i9 



47r 
+ sin"^ 1? 



1 + cos(q ■ d)e 



1)?) 



}■ (40) 



The fringe visibility in this case can approach 100% in the 
i9 = 7r/2 plane if the Debye- Waller factor is close to 1. 
The cross section for detection of light with polarization 

is 



= 2 



dr^out 

which is totally isotropic and shows no fringes. 

D. Which-path interpretation 

The presence of interference fringes in the tt case and their 
absence in the a case has a simple explanation in terms of the 
possibility, in principle, of determining which of the two ions 
scattered the photon. Consider the sequence of transitions in 
Fig. |^(a), representing the tt case. Each box represents the 
combined state of the two ions. Ion 1 is represented by the 
diagram on the left side of a box and ion 2 by that on the 
right. The ordering of energy levels is the same as in Fig. ^. 
For simplicity, we neglect the translational degrees of freedom, 
which lead to the appearance of the Debye- Waller factor in 
Eq. (P9). The system begins in the state 



I*.) = l('Si/2, +1/2)1(^51/2, -1/2)2). 



(42) 



One ion or the other absorbs a photon from the laser beam and 
undergoes a tt transition to the excited state. That ion emits 
a photon and undergoes a tt transition back to the ground 
state. The two paths, corresponding to either ion 1 or ion 2 
scattering the photon, lead to the same final state. Therefore, 
the amplitudes for these two paths must be added, and this 
leads to interference. Since the final states of the ion are the 
same as the initial states, it is not possible to determine which 
of the ions scattered the photon by examining their states. 

Now consider the sequence of transitions in Fig. ^(b), rep- 
resenting the a case. As in the previous case, one ion or the 
other absorbs a photon and undergoes a tt transition to the 
excited state. However, in this case, that ion undergoes a a 
transition when it emits a photon and changes its mj quan- 
tum number. The final states differ, depending on which of 
the ions scattered the photon. Hence, there is no interference 
between the two paths. It would be possible to tell which 
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ion scattered the photon by examining the states of the ions 
before and after the scattering. 

The preceding analysis is valid only in the limit of low laser 
intensity, so that the probability of both ions being excited 
at the same time is negligible, and stimulated emission can 
be neglected. It is not necessary that the two ions be in the 
same quantum state for interference to occur, only that the 
final combined states for the two paths be indistinguishable. 
For definiteness, a particular initial state [Eq. (^^] was cho- 
sen. For each of the 3 other possible initial states, there is 
a process like Fig. ^(a) in which the ions scatter a photon 
and return to their original states and one like Fig. ^(b) in 
which one of them scatters a photon and changes its state. 
Processes of the former type lead to interference; those of the 
latter type do not. 



VI. COMPARISON WITH EXPERIMENT 

Figure |^ shows an image of the fringes observed for the vr 
case, in which ein was perpendicular to the X-Z plane and 
the detector was sensitive only to light polarized parallel to 
ein. The dark spots are due to stray reflections of the laser 
beams. When ein was rotated by 90° without changing the 
polarizer in front of the detector (a case), the image showed 
no fringes. The image data from a single ion, which shows no 
interference fringes, were used to correct the data of Fig. |^ 
for a slowly spatially varying detection efficiency. The data 
within the rectangle in Fig. H were summed along the vertical 
direction and divided by the detection-efficiency function. 

The normalized data points are shown in Fig. |^ together 
with a least-squares fit. In this fit, as in Ref. the temper- 
atures of the stretch and tilt modes were assumed to have the 
ratio expected from theory , 

TtiTt = {1 + [3cos^(e)]-^}/{l + [3sin^(e)]-^}, (43) 

and both temperatures were allowed to vary together in the 
fit. The fringe visibility in the vicinity of the X-Z plane 
is insensitive to the temperature of the Y motion, which is 
cooled indirectly by coupling to the other modes. The mean 
ion separation was calculated from knowledge of the trap pa- 
rameters. The dependence of Eq. on the out-of-plane 
angle "I> is small, and <1> was set to in the fit. The fitted 
value of T"^^ was 1.08 ± 0.12 mK, or 0.92 ± 0.10 times the 
Doppler-cooling limit. The fringe visibility, extrapolated to 
(/> = 0, would be 100% if it followed Eq. (|2|). The fitted value 
for this parameter was (71 ± 4)%. The errors represent the 
standard deviations estimated from the fit. The maximum 
observed visibility, at the minimum value of 4> in Fig. ^ is 
approximately 60%. 

There are several likely causes of the difference between the 
observed and predicted values of the fringe visibility. First, 
the theory was derived for the limit of low intensity. The 
saturation parameter was measured to be s = 0.078 ± 0.025 
(see Appendix). By itself, this would reduce the maximum 
visibility to (1 -I- s)~^ « 93%, because the spectrum of the 
resonance fluorescence in this polarization contains an inco- 
herent part Other likely causes of reduced visibility are 
unequal laser intensities at the two ions, imperfect polarizers, 
stray background light, and quantum jumps of one of the ions 



to a metastable state, leaving only one ion fluorescing. Each 
of these effects might reduce the visibility by a few percent. 



VII. DISCUSSION 

The fact that the resonance fluorescence from a two-level 
atom illuminated by weak, monochromatic light is coherent 
with the applied fleld was noted by Heitler . The spectrum 
of the resonance fluorescence for arbitrary applied intensities 
was calculated by MoUow |^ . In the limit of low applied in- 
tensity, the spectrum is monochromatic and coherent with the 
applied field (a S function). At higher intensities, the coherent 
component decreases in amplitude, and a component not co- 
herent with the applied field and having a width equal to the 
natural linewidth appears. At very high intensities, the co- 
herent component continues to decrease in amplitude, and the 
incoherent component splits into three separate Lorentzians. 
The existence of a coherent component in the resonance flu- 
orescence of a single ion was confirmed directly by Hoffges et 
al. by a heterodyne measurement |p3| . 

Classically, we would expect the resonance fluorescence 
from two two-level atoms at fixed positions, excited by the 
same monochromatic field, to generate interference fringes 
having 100% visibility in the limit of low applied intensity, 
since the radiated fields are coherent with each other. At 
higher applied intensities, the visibility should decrease, since 
more of the resonance fluorescence intensity belongs to the in- 
coherent component. Quantum treatments for two two- level 
atoms have been given by Richter and by Kochan et at 
p5[ , who predict a visibility equal to (1 -I- s)~^, where s is 
the saturation parameter defined in Ref. |2^. This is just the 
ratio of the intensity of the coherent component to the total 
resonance fiuorescence intensity for a single atom. 

Polder and Schuurmans calculated the spectrum of the 
resonance fiuorescence ofaJ = l/2toJ = l/2 transition for 
a single atom. The spectrum of the light having polarization 

is like that for a two-level atom. Hence, interference fringes 
would be expected in the e.n--polarized resonance fluorescence 
from two such atoms for low applied intensity. The spectrum 
of the light having polarization does not contain a S func- 
tion. In the limit of low applied intensity, it is a Lorentzian 
having a width approximately equal to the photon scattering 
rate, which can be much less than the natural linewidth. Even 
for applied intensities approaching s = 1 , the coherence length 
is on the order of c/7, where 7 is the spontaneous decay rate 
of the excited state. For the Hg"*" 6p^Pi/2 level, this is about 
70 cm. For interference fringes to exist, the radiation from the 
two atoms must be mutually coherent. Whether or not fringes 
should exist in the eo--polarized light from two atoms is not 
immediately obvious from a classical analysis. However, the 
perturbative quantum treatment of Sec. V predicts that there 
should be no interference, since there is only one probability 
amplitude connecting the initial and final states. The absence 
of interference in this case is fundamentally a quantum effect, 
though one having more to do with the quantum nature of 
the atom and the existence of degenerate, orthogonal ground 
states, than with the quantum nature of the electromagnetic 
field. Precisely the same point was made by Scully and Driihl 
when they showed that interference fringes are not present in 



7 



the Raman radiation emitted by two three-level atoms having 
a A configuration l^l. 

Wong et al. (hsT calculated the interference of resonance 
fiuorescence from two four-level atoms having a level struc- 
ture like that of ^^^Hg"*". Their analytic calculations are for a 
simpler geometry than the one actually used by Eichmann et 
al. and ignore the motion of the ions. They do, however, 
include the effect of the decrease in visibility due to the inco- 
herent component of the resonance fiuorescence, which is not 
included in the perturbative calculation of Sec. V. The ana- 
lytic calculations of Wong et al. and the present calculations 
agree in the limits in which they are both valid, that is, for low 
applied intensities and for no ion motion. In particular, Wong 
et al. show that the fringe visibility can approach 100% at 
low applied intensities, with polarization-selective detection. 
Wong et al. also made Monte Carlo wavefunction simulations, 
in which the motion of the ions was included classically, and 
observed a decrease in visibility due to this effect. 

Huang et al. calculated the effect of thermal motion 
on the interference fringe visibility for two two-level atoms, 
each trapped in a separate harmonic well. They obtained an 
expression equivalent to Eq. (1) of Eichmann et al. |Q] for this 
model. However, the treatment of Eichmann et al., the de- 
tails of which are given in the present article, is more useful 
for the analysis of the experiment of Ref. [Q], since it deals 
explicitly with the actual normal mode structure of the two 
trapped ions. 

Brewer has published a theory of interference in the light 
scattered from two four- level atoms jll]. One prediction of 
this theory is that the fringe visibility cannot exceed 50%, 
even with polarization-selective detection. This contradicts 
the experimental results of Sec. VI shown in Fig. ^. While 
the maximum visibility is about 60%, only slightly exceed- 
ing 50%, no background has been subtracted from the data, 
and there are several known sources of decreased visibility, in- 
cluding thermal motion of the ions, the incoherent component 
of the resonance fluorescence, and stray scattered light. The 
data were normalized by division by a slowly varying detec- 
tion sensitivity function, a process that cannot enhance the 
visibility. 

The basic flaw in Brewer's argument can be seen in Eq. (2) 
of Ref. where he lists the basis states for the two-atom 
system. The states |5)-|8) are the four states in which both 
atoms are in the ground electronic state. The states |1)-|4) 
are linear combinations of states in which one atom is in the 
ground state and one is in the excited state. However, most 
of the possible states of this type are missing, apparently be- 
cause of a false assumption that the allowed states must have 
a particular kind of exchange symmetry. For example, the 
intermediate superposition state shown in Fig. ^a) is, in his 
notation. 



(|ci62)7r„„,i -I- |aid2)7r„„,2), 



(44) 



and is not contained in the list. The neglect of these basis 
states leads to the neglect of processes like that of Fig. ^a), 
in which the two atoms are initially in different mj states. 
Thus, he reaches the false conclusion that the two atoms must 
initially be in the same mj state in order for interference to 
occur. Since he misses half of the processes that lead to inter- 



ference, he predicts a maximum visibility, with polarization- 
sensitive detection, of 50% rather than 100%. 

We conclude with some remarks regarding the principle of 
complementarity. Wave and particle properties of light are 
complementary and hence cannot be observed at the same 
time. If it is possible to determine which atom scattered 
the photon, the interference fringes must vanish. Feynman's 
thought experiments, in which various methods of determin- 
ing the path of an electron through a two-slit Young's inter- 
ferometer lead to the destruction of interference fringes due 
to a random momentum transfer are often quoted [Ref. 
pp. 1-6-1-11]. However, in Ch. 3 of the same textbook, Feyn- 
man emphasizes the seemingly more fundamental viewpoint 
that interference is present only if there exist different indis- 
tinguishable ways to go from a given initial state to the same 
final state. His example of the scattering of neutrons from a 
crystal is very similar to the experiment of Eichmann et al. If 
the nuclei of the atoms in the crystal have a nonzero spin, the 
angular distribution of scattered neutrons is the sum of a fea- 
tureless background and some sharp diffraction peaks. The 
sharp diffraction peaks are associated with neutrons which do 
not change their spin orientations in the scattering. The fea- 
tureless background is associated with neutrons whose spins 
change their orientations in the scattering. In this case, there 
must also be a change in the spin orientation of one of the 
nuclei in the crystal. It would be possible in principle to de- 
termine the nucleus which scattered the neutron, so there is 
no interference. The position-momentum indeterminacy re- 
lations play no essential role in the presence or absence of 
interference in this case, just as in the experiment of Eich- 
mann et al. 
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APPENDIX: CALIBRATION OF THE 
SATURATION PARAMETER 

For the case where an electric-dipole transition between a 
ground state and a ^Pi/2 excited state is excited by 
linearly polarized light, we define the saturation parameter 
s similarly to the way in which it is defined for a two-level 
system |26|. The magnetic field is assumed to be small, and 
the quantization axis for the ion is along the electric field. We 
define 



(7/2)2 



(Al) 



where Q.i = 6"^/^|£:o(^S'i/2 |pPi/2)|?i~^ is the Rabi fre- 
quency, and the other terms have been defined previously. 
In order for the perturbative analysis of Sec. V to be valid, 
we must have s ^ 1. In the case of Hg+, the ^P\/2 state 
has a small (approximately 10~^ ) probability of decaying 
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to the metastable ^-Da/g state, which decays either directly 
to the ground state or to the metastable ^1)5/2 state, which 
decays to the ground state. The 194.2 nm fluorescence inten- 
sity from a single ion is bistable, since it has a steady level 
when the ion is cycling between the ^'id '^P\/2 states 

and vanishes when the ion drops to a metastable state. The 
fractional population of the state, summed over both 

mj values, is s/[2(l + s)] while the ion is cycling between the 
Si/2 and Pi/2 states. The quantum jump statistics have 
been discussed in several previous articles [ p^p^ , p8t . For a 
single ion, we define Pon to be the fraction of the time that 
the ion is cycling between the ^Si/2 and ^P\/2 states, and 
Poff=(l — Pon) to be the fraction of the time that it spends 
in either of the metastable states. It can be shown, from 
the steady-state solutions of the differential equations for the 
populations [Eqs. (2a)-(2c) of Ref. |2^]], that s is related to 
the ratio Pos/Pon according to 

1 S _ 7l72(Poff/pon) ^ PoS ,.„s 
2{1 + S)- 73(72 +/271) ^ Pon' ^ ^ 

where the parameters 71, 72, 73, and /2 have been measured 
[ p^ , and the uncertainty in the coefficient (0.36) is about 30%, 
due mostly to the uncertainty in 73. 

For two ions, the fluorescence will be tristable, since 0, 1, 
or 2 ions may be in a metastable state. During an interference 
fringe measurement, the number of photons detected in each 
successive period of a few milliseconds was recorded. Fig- 
ure ^ shows a plot of the probability distribution of the 5 ms 
photon counts during the measurement of Fig. ^ The three 
peaks correspond, from left to right, to 2, 1, or ions being in 
a metastable state. The leftmost peak corresponds to the sig- 
nal from stray background light, since there is no fluorescence 
from the ions. The curve is a least-squares fit to a sum of three 
Gaussians. The areas under the peaks should be in the ratio 
Poff :2poffPon:Pon- The ratios of the areas obtained from the 
fit are 0.011:0.160:0.828, so Poff /Pon=0.10 ± 0.01, and, from 
Eq. ( A2), s=0.078±0.025, so the perturbative analysis should 



be a good approximation. 

During this measurement period, the interference fringe de- 
tection was gated off for 5 ms if the number of photons de- 
tected in the previous 5 ms was less than 80. This helped 
to prevent loss of the fringe visibility due to background 
from single-ion fluorescence, which would have no interference 
fringes. 
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FIG. 1. Geometry of the Young's interference experiment, 
projected onto the X-Z plane. The equilibrium positions of 
the two ions, represented by the filled circles, lie along the 
Z axis. The wavovcctor kin of the incoming photon is in 
the X-Z plane, making an angle O with the Z axis. The Y 
axis is out of the plane of the figure. The projection of the 
wavevector kout onto the X-Z plane makes an angle </> with 
kin. The angle that kout deviates from the X-Z plane in the 
+y direction is $ (not shown). The polarization vectors of 
the incoming and outgoing photons are ein and €out. 
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FIG. 2. Zeeman sublevels involved in the 194 nm, 6s^iSi/2 
to 6p^Pi/2 transition of ^®*Hg^. The allowed tt and a tran- 
sitions are labeled. The Zeeman splitting of the levels is ex- 
aggerated. 



FIG. 3. Coordinate system for description of the direction 
and polarization of the outgoing photon. The z axis is parallel 



to Cin, and the x axis is parallel to kin 
vector e„ lies in the plane containing ein 
is perpendicular to that plane. 



The polarization 
and kout, while en- 
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(a) 



(b) 



FIG. 4. Each box represents the combined state of the two ions. The ordering of energy levels is the same as in Fig. ^. 
In (a) (the tt case), one ion or the other undergoes a tt transition from the ground to the excited state. That ion undergoes 
a TT transition back to the ground state. The two paths lead to the same final state of the two ions. Hence, the probability 
amplitudes must be added, and interference is possible. In (b) (the a case), one ion or the other undergoes a tt transition to 
the excited state, but the excited ion undergoes a a transition to the ground state. The two paths lead to different final states 
of the two ions. Hence, there is no possibility of interference. In order for interference to occur, it is not necessary that the 
initial states of the two ions be the same, only that the final combined states for the two paths be the same. 
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FIG. 5. Experimental fringe data for the case in which the 
the detected light is polarized in the same direction as the in- 
coming light (tt case). The ion separation d = 4.17 ^m. The 
angle 4> (the deviation from the forward scattering direction) 
increases to the right. The decrease in visibility with increas- 
ing <j) is due to thermal motion of the ions. The dark spots are 
due to stray reflections of the laser beams. The data within 
the rectangle were summed along the vertical direction and 
least-squares fitted. 



I I 

20 25 30 35 40 45 50 

FIG. 6. Experimental fringe data (dots) from the image of 
Fig. P and a least-squares fit (line) to the sum of the theoreti- 
cal intensity [Eq. (Ea)] and a constant background. The fitted 
temperature is approximately equal to the Doppler-cooling 
limit. 
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FIG. 7. Plot of the probability distribution of the fluores- 
cence intensity for two ions, used to determine the saturation 
parameter s. The horizontal axis corresponds to the number 
of photons counted in a 5 ms interval. The vertical axis cor- 
responds to the number of 5 ms intervals in which a given 
number of photons was counted. This was measured simulta- 
neously with the interference fringes shown in Fig. ^ The curve 
is a least-squares fit to a sum of three Gaussians. The areas 
under the Gaussians, from left to right, are proportional to 
the probabilities that 2, 1, or of the ions are in a metastable 
state. Higher values of s correspond to higher populations in 
the metastable states. 



12 



